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SYSTEME D' INTELLIGENCE ARTIFICIELLE POUR LA 
CLASSIFICATION D' EVENEMENTS , D' OBJETS OU DE SITUATIONS 
A PARTIR DE SIGNAUX ET DE PARAMETRES DISCRIMINANTS 

ISSUS DE MODELES 

5 

DESCRIPTION 

DOMAINE TECHNIQUE 

La presente invention concerne un systeme 
10 d f intelligence artificielle pour la classification 
d'eyenements, d'objets ou de situations a partir de 
signaux et de parametres discriminants issus de 
modeles.. 

L' invention s' applique notamment a la 
15 classification d' evenernents sismiques. Une telle 
classification sera consideree a titre d'exemple non 
limitatif dans la suite de la description. 

ETAT DE LA TECHNIQUE ANTERIEURE 

20 La classification automatique d / evenernents sismiques 

La classification automatique d' evenernents 
sismiques est un probleme relativement recent, puisque 
le probleme n' a veritablement ete aborde que dans les 
annees 1980. Ces travaux sont maj oritairement orientes 

25 vers la recherche de parametres discriminants (c'est-a- 
dire susceptibles de perrnettre la classification) dans 
les signaux sismiques. De nombreuses caracteristiques 
potentielles sont proposees en vue d'une future 
classification automatique. A partir de 1990, on 

30 commence a trouver dans les articles publies des. 
tentatives de classification automatique, que ce soit 
par des techniques neuronales ou par des systemes a 
base de regies. Ces travaux cherchent a separer les 
seismes d' origine naturelle des explosions. Aucun de 
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ces articles ne traite de la discrimination 
d' ef f ondrements miniers coups de terrain ") . 

Du fait de la complexity du probleme, ces 
articles mettent clairement en evidence la necessite de 
5 mettre au point des systemes automatiques capables 
d' apprentissage. Les methodes neuronales ont done 
souvent ete proposees pour r6aliser la discrimination 
automatique d'evenements sismiques, mais avec des 
limitations qui sont analysees plus loin. Les modeles 
10 les plus couramment proposes sont les perceptrons 
multi-couches a connexions completes entre couches 
successives . 

Tous ces articles cherchent a determiner 
l'origine du seisme a partir de caracteristiques 
15 extraites uniquement des signaux. Les donnees de plus 
haut niveau (date, heure, latitude, longitude, 
magnitude...) ne sont jamais employees pour la 
classification. Or, les sismologues connaissent la 
difficulty de discriminer les signaux sismiques 
20 uniquement par des traitements de bas niveau. 

Les travaux de Baumgardt et de ses 
collaborateurs decrits dans le document reference [1] 
sont sans doute ceux qui ont permis les plus grandes 
avancees dans la recherche de parametres 
25 discriminants. 

Les variations du cepstre, le cepstre d' un 
signal x etant la transformer de Fourier inverse du 
logarithme de la transformee de Fourier de x, sont 
souvent utilisees. On peut ainsi montrer que le cepstre 
30 permet de visualiser le phenomene de micro-retards 
presents dans les signaux de tirs caracterises par une 
plus large variance. Le document reference [2] note 
aussi cette propriete en signalant neanmoins que 
1' absence de cette caracteristique ne permet aucune 
35 deduction quant a la classe de l'evenement. 
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Les rapports des amplitudes des differents 
types d'ondes peuvent aussi servir de discriminants. Le 
document reference [3] etudie toute une serie de 
rapports d' amplitudes (P n /L g , P g /L g/ L g /R g ) . Ces rapports 
5 sont annonces comme pouvant permettre une bonne 
discrimination. 

Les memes auteurs introduisent aussi les 
rapports des densites spectrales de puissance des 
differents types d'ondes detectes. Comme pour les 

10 rapports d' amplitudes , ces discriminants sont utilises 
par toutes les etudes qui recherchent des discriminants 
dans les signaux sismiques. Pour caracteriser les 
explosions, on a egalement utilise les rapports des 
densites spectrales de puissance d'un type d'onde, ici 

15 S, sur des bandes de frequences differentes, c'est-a- 
dire le rapport de la densite spectrale de puissance de 
S dans la gamme 1-2 Hz a la densite spectrale de 
puissance de cette meme phase dans la bande 7-20 Hz. Le 
rapport entre des densites spectrales de puissance de 

20 l'onde S avant et apres 10 Hz est egalement presente 
comme un bon separateur des explosions des seismes. 

Le document reference [4] remarque que le 
temps de propagation des signaux provenant d'une mine a 
un temps t Sg -t Pg constant pour une station 

25 d' enregistrement donnee. Ce temps de propagation est 
presente comme une caracteristique potentielle pour une 
mine, restant cependant moins fiable que les 
precedentes caracteristiques . 

Le document reference [5] propose d' utiliser 

30 la presence de l'onde de surface des seismes pour les 
discriminer des explosions nucleaires a des distances 
regionales. La caracterisation de la presence d'une 
onde de surface est effectuee indirectement en 
comparant les magnitudes m b et M s . Pour deux evenements 

35 sismiques de meme magnitude m b , la magnitude de l f onde 
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de surface M s es.t en general plus 61evee dans le cas 
d'un seisme a cause de la presence de l'onde de surface 
que dans le cas d'une explosion. En effet, cette onde 
de Rayleigh de croute intervient dans le calcul de la 
5 magnitude M s et sa presence est subordonn6e au 
phenomene de cisaillement absent dans le cas 
d' explosions nucleaires. La representation de la 
difference (m b - M s ) en fonction de m b permet de 
verifier cette hypothese. Neanmoins, le calcul de la 

10 magnitude M s depend de la periodicite du signal 
enregistre et n' est pas possible en toute rigueur pour 
des evenements regionaux. Au contraire, la presence 
d'une onde de surface, correspondant a une onde de 
Rayleigh sedimentaire dans des signaux sismiques 

15 proches, caracterise des evenements de nature 
artif icielle . Une methode de detection de ce second 
type d' onde de surface consiste a rechercher 
directement sa presence dans le spectrogramme du 
signal, sa frequence etant connue (entre 0,5 et 1,5 Hz) 

20 et son heure d' arrivee supposee pouvant etre calculee a 
partir de sa vitesse moyenne de propagation et de la 
distance separant 1' epicentre de la station 
d' enregistrement . 

25 Les systemes decrits dans les documents de 

l'art connu ne sont pas operationnellement credibles 
pour plusieurs raisons : 

• Les etudes menees par des geophysiciens, le 
plus souvent riches et detaillees en matiere de 

30 proposition de parametres discriminants, ne proposent 
pas de methode fiable pour 1' exploitation automatique 
de ces parametres. 

• Les etudes menees par des 
inf ormaticiens proposent des systemes qui ne tirent pas 
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assez parti de la complementarite des donnees et des 
connaissances geophysiques . 

La plupart des etudes de 1' art connu 
utilisent des bases de donnees d' ev^nements sismiques 
5 de taille extremement reduite, ce qui a pour 
consequence de ne pas permettre un apprentissage 
statistique correct. La classification se fait, le plus 
souvent, sur des bases ayant moins d' une centaine 
d'evenements, comme decrit dans les documents 

10 references [2] et [6]. L'une des bases les plus 
importantes trouvee dans les documents de l'art connu 
est composee de 312 evenements seulement, comme decrit 
dans le document [4] . La consequence directe en est que 
les marges d'erreurs sur les resultats presentes sont 

15 tres elevees, ce qui ne permet pas d'accorder beaucoup 
de confiance k ces resultats. 

La repartition geographique des exemples de 
la base de donnees est un element tres important. La 
plupart des bases regroupent des evenements qui ont eu 

20 lieu dans des regions de taille restreinte (quelques 
dizaines de kilometres de cote), ou les proprietes 
geologiques des sous-sols sont peu diversif iees . La 
recherche de discriminants generaux est done biaisee, 
les discriminants n'etant efficaces que pour une region 

25 donnee. 

De plus, comme decrit dans le document 
reference [1], les evenements des deux classes a 
discriminer peuvent provenir de deux regions 
geographiques clairement distinctes, parfois eloignees 

30 l'une de 1' autre de plusieurs centaines de kilometres. 
II est alors impossible de savoir dans quelle mesure ce 
n'est pas la " coloration " des signaux par les couches 
geologiques traversees qui permet la discrimination, 
plus que les signaux eux-memes. Or les sismologues 

35 savent que cette coloration est loin d'etre negligeable 
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et que 1' information de localisation est tres 
importante. 

Le plus souvent, un nombre tres limits de 
stations d' enregistrement est utilise. Les signaux sont 
5 enregistres par deux ou trois stations au maximum, mais 
le plus souvent on se contente d'une unique station. 
L'6venement sismique est alors represents par un seul 
signal, ce qui reduit considerablement 1' information 
utilisable. 

10 Enfin, les evenements integres dans la base 

de donnees sont tres generalement selectionnes selon un 
critere prealablement defini : magnitude superieure a 
un seuil, rapport signal/bruit superieur a un seuil, 
comme decrit dans le document reference [2]. Mais cette 

15 selection biaise evidemment totalement les resultats. 

Si les etudes de l'art connu ont permis de 
recenser une large gamme de discriminants 
potentiellement utiles a la classification, il s'est 
av6re tres difficile de trouver des discriminants 

20 globaux efficaces en raison du nombre eleve de types de 
tirs et de seismes distincts. 

La plupart du temps, les classifieurs 
proposes se revelent inefficaces parce que trop simples 
(separateurs lineaires) ou impossibles a regler du fait 

25 de leur complexite. lis imposent un enorme travail de 
pretraitement sur les donnees, ce qui rend les syst ernes 
proposes non generalisables . 

Exemple de surveillance sismique 

30 Le laboratoire de detection et de geophysique 

(LDG) du CEA surveille en continu 1' activite sismique 
de la terre depuis 1962. Lorsqu'un evenement sismique 
survient en un endroit du globe, il est enregistre en 
France par un reseau de quarante-deux sismometres 

35 verticaux localises sur le territoire metropolitain, 
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comme illustre sur la figure 1, les stations CP etant 
des stations courtes periodes, et les stations LP des 
stations longues periodes. Une description detaill§e du 
reseau de sismometres et de la propagation des ondes 
5 sismiques en France est donnee dans le document 
reference [7]. 

Ce reseau, qui utilisait depuis sa creation 
la transmission par voies hertziennes, est passe 
recemment a une transmission numerique par satellite. 

10 Un filtrage et un gain adaptes du signal permettent de 
detecter les seismes proches ou au contraire les 
seismes lointains de plus longue periode appeles 
teleseismes. Le reglage des parametres de filtrage et 
du gain doit permettre de trouver un compromis entre la 

15 detection des evenements sismiques de magnitude 
relativement faible et le bruit de fond. 

La figure 2 montre des signaux enregistres 
par les sismometres du laboratoire LDG situes entre 84 
et 14 6 kilometres (soit SBF : 84 km ; PGF : 110 km ; 

20 FRF : 127 km ; LMR : 136 km et LRG : 148 km) de 
l'epicentre estime d'un seisme de magnitude 1,9 
localise a 10 kilometres au sud d'Imperia en Italie, le 
9 mai 1996 (heure : 1 heure O minute 59 secondes ; 
latitude : 43.34 ; longitude : 8.19 ; magnitude : 1,9). 

25 Sur chaque signal sont repertoriees 

differentes phases sismiques qui serviront lors de 
1' analyse detaillee de l'evenement. 

Chaque annee, environ 9000 evenements 
sismiques proches sont ainsi detectes, dont 800 a 1 200 

30 seismes naturels. 

Les sismologues du laboratoire LDG 
depouillent et analysent quotidiennement les donnees 
enregistrees par les stations du reseau frangais. lis 
publient un bulletin hebdomadaire contenant 1' ensemble 

35 des seismes naturels survenus en France ou dans les 
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regions limitrophes . Un bulletin similaire est publie 
pour les teleseismes. Le tableau I, donne en fin de 
description, est un extrait du bulletin relatif £ la 
periode du 9 au 15 septembre 1998. Dans ce tableau on a 
5 les abreviations suivantes : 

HEURE OR : heure d' origine (TU) ; 

LAT : latitude de 1' epicentre (deg) ; 

LON : longitude de 1' epicentre (deg) ; 

PRO : profondeur de 1' epicentre (km) ; 

10 ML : magnitude locale ; 

RMS : residu quadrat ique mo yen (s) . 

Ce tableau comporte done l f ensemble des s^ismes proches 
detectes durant cette semaine et les caracteristiques 
de chaque evenement : date et heure origine, position 

15 et profondeur de 1' epicentre, magnitude, r6sidu 
quadratique moyen de localisation et region de 
localisation deduite. 

En raison du fort desequilibre entre le 
nombre d' evenements artificiels et celui d' evenements 

20 naturels (en France, les evenements sismiques d' origine 
artificielle sont dix fois plus frequents que les 
evenements sismiques naturels) , seuls les evenements 
supposes etre des seismes ou les evenements de classe 
indeterminee sont extraits du bruit de fond par les 

25 sismologues pour etre ensuite analyses plus precisement 
par un logiciel de localisation. Les autres signaux 
(evenements artificiels pour la plupart) sont archives 
pendant six mois. 

La procedure d' exploitation des signaux 

30 consiste en une localisation deja automatisee, suivie 
d'une phase de caracterisation (determination de 
1' evenement a 1' origine des signaux). 

L' exploitation des signaux est effectuee a 
l'aide du systeme global illustre sur la figure 3. 

35 
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1/ invention a pour, objectif de pallier les 
inconvenients des systemes de l'art connu en proposant 
un nouveau systeme d' intelligence artificielle pour la 
classification d' ev6nements, d'objets ou de situations 
5 a partir de signaux et de parametres discriminants 
issus de modeles. 

EKPOSE DE 1/ INVENTION 

La presente invention concerne un systeme 

10 d' intelligence artificielle pour la classification 
d'evenements, d'objets ou de situations a partir de 
signaux et de parametres discriminants issus de 
modeles , caracterise en ce qu'il comprend au moins une 
branche de traitement comprenant un systeme expert flou 

15 (expert SEF) prenant une decision a partir de 
proprietes de haut niveau et de parametres 
discriminants de plus bas niveau extraits des signaux 
par des procedures de type traitement de signal, et 
capable d'expliquer sa decision a 1' utilisateur par 

20 1' intermediaire de regies triees par ordre 
d' applicability . 

Dans ce systeme expert flou on effectue une 
descente de gradient sur les parametres : 
x = y/a 

25 • s = ln(2a 2 ) 

• r = ln(p) 

• d 

avec : 

y : position des ensembles flous des 
30 premisses 

a : largeur des ensembles flous des 
premisses 

p : poids des regies 

d : degre d' activation de chaque classe 
35 pour chaque regie 
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Avantageusement le systerne de 1' invention est 
un systerne multi-experts constitue d'au moins deux 
branches de traitement independantes, se conf igurant 
automatiquement par apprentissage statistique sur des 
5 bases de donnees, ayant des proprietes particulieres, 
et fusionnees par un systerne decisionnel de haut 
- niveau- Avantageusement une branche comprend un 
classifieur neuro-flou (expert CNF) prenant ses 
decisions a partir de proprietes de haut niveau et de 
10 parametres discriminants de plus bas niveau extraits 
des signaux par des procedures de type traitement de 
signal. Avantageusement une autre branche comprend un 
reseau neuronal a connexions locales et poids partages 
(expert TDNN) constitues de bancs de filtres adaptatifs 
15 non lineaires, extrayant lui-meme des informations 
discriminantes de representation temps-frequence des 
signaux correspondants a l'ev6nement. 

L' invention peut etre utilisee dans 
dif f erents domaines d' application et notamment en : 

20 

# Suzrrelllance d ' erenements geophy&lqaes 

Le systerne s' applique alors a 1' analyse de 
tout evenement geophysique observable par des signaux 
recueillis par des stations : 
25 - signaux sismiques ; 

- infrasons ; 

- ondes hydro-acoustiques . 

Ces evenements peuvent etre des evenements 
proches (dits regionaux) ou lointains (par exemple des 
30 teleseismes) . 

la fonction a assurer peut etre : 

- un filtrage pour eliminer les evenements 
non pertinents pour des traitements en aval ; 

- une detection d' evenements particuliers ; 
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- une classification exhaustive en un 
ensemble de groupes d' evenements de meme nature. 

# Surveillance et con t role industrlels 

5 Le systeme s' applique egalement a 1' analyse 

d'objets ou de processus industriels, pour peu que l'on 
dispose de signaux ou d' images recueillis par des 
capteurs. En voici quelques exemples : 

■ Controle de qualite d'objets ou de produits 
10 manufactures : l'objectif est de controler la forme 

et/ou la position d'objets, detecter et caracteriser 
des defauts. Les experts CNF et SEF utilisent des 
mesures realisees par traitement d' images. L' expert 
TDNN utilise une ou plusieurs images de la piece. 

15 ■ Maintenance predictive d' equipements : 

l'objectif est de prevoir une future defaillance de 
machines, d' ordinateurs, d' equipements electroniques, 
de capteurs afin de donner l'alerte et de pouvoir 
mettre en oeuvre une procedure de correction avant la 

20 panne. Les experts CNF et SEF utilisent des mesures de 
coefficients de haut niveau, des correlations. L' expert 
TDNN utilise des signaux. 

■ Surveillance de processus complexe : 
l'objectif est de surveiller le bon f onctionnement 

25 d'une chaine de production. Les experts CNF et 
SEF utilisent des mesures de coefficients de haut 
niveau. L' expert TDNN utilise des mesures de 
coefficients de bas niveau. 



30 Dans le domaine geophysique, les proprietes 

de haut niveau peuvent etre la localisation, la 
magnitude, l'heure et la date. Le systeme de 
1' invention permet de realiser la classification 
automatique d' evenements sismiques en trois ensembles : 

35 - seismes naturels ; 
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- explosions (tirs de mines et essais 
militaires) ; 

- coups de terrain (ef f ondrement de couches 

minieres) . 

5 Le systeme s' integre alors dans une chaine de 

traitements automatiques pour op§rer une fonction de 
filtrage des evenements sismiques. Ses caracteristiques 
principales sont : 

- une fiabilite maximale : le systeme est 
10 capable de prendre des decisions meme avec des donnees 

degradees ou imprecises, voire en 1' absence de 
certaines informations ; 

- l'acces a 1' explication des decisions afin 
de lever un doute eventuel sur une decision. 

15 

BREVE DESCRIPTION DES DESSINS 

La figure 1 illustre le reseau de sismometres 
du laboratoire de detection et de geophysique (LDG) du 
CEA en 1998 . 

20 La figure 2 illustre des exemples de signaux 

sismiques enregistres par le reseau de la figure 1. 

La figure 3 illustre un systeme global 
d' exploitation de signaux geophysique de l'art connu. 

La figure 4 illustre le schema de principe du 
25 systeme multi-experts de discrimination d' evenements 
sismiques de 1' invention. 

La figure 5 illustre le schema general 
d' apprentissage du systeme de 1' invention a partir 
d' exemples . 

30 La figure 6 illustre un neurone artificiel. 

La figure 7 illustre un reseau de neurones 
artif iciels . 

La figure 8 illustre un classifieur neuro- 

flou. 
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La figure 9 illustre un mecanisme 
d' activation de cellules de codage. 

La figure 10 illustre des exemples 
deactivation des cellules de codage. 
5 La figure 11 illustre un systeme expert flou 

pour la discrimination d'evenements sismiques. 

Les figures 12A a 12D illustrent des 
pretraitements successifs appliques a un signal 
sismique. 

10 La figure 13 est une architecture d'un reseau 

de neurones a connexions locales et poids partages. 

La figure 14 illustre la representation des 
epicentres des evenements sismiques de 1962 a 1996. 

15 EXPOSE DETAILLE DE MODES DE REALISATION 
Description generale du systeme 

Le systeme de 1' invention comporte au moins 
une branche de traitement contenant un systeme expert 
flou. Lorsqu'il comprend plusieurs branches de 

20 traitements independantes, on parle de systeme multi- 
experts. 

• Prise de decision multi-experts 

Le principe de la prise de decision multi- 
25 experts, qui est done 1'un des modes de realisation de 
1' invention, est 1' exploitation de la synergie entre 
plusieurs branches de traitement complementaires . Cette 
complementarite reside dans : 

A. Les performances generales : une branche 
30 est plutot generaliste (performances assez 

bonnes dans la plupart des cas) , une autre 
est plutot specialiste (tres bonnes 
performances sur certains cas difficiles, 
taux d'erreurs plus eleve sur les cas ne 
35 relevant pas de sa competence) . 
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B. Les performances selon le cas traite : une 
branche peut etre plus apte qu' une autre a 
traiter tel cas particulier. 

C. La nature des entrees (signaux ou donnees 
5 de haut niveau) . 

D. La nature des sorties (simple donnee de la 
classe, estimation de la certitude de la 
decision, explication formelle de la 
decision) . 

10 La figure 4 illustre le schema de principe du 

systeme multi-experts de discrimination d' evenements 
sismiques de 1' invention. Ce systeme est constitue de 
plusieurs branches de traitement independantes, chacune 
ayant des proprietes particulieres , fusionnees par un 

15 systeme decisionnel de haut niveau. 

Ces branches sont : 

- un classifieur neuro-flou, dit CNF, prenant 
ses decisions a partir de proprietes de haut niveau des 
evenements (par exemple pour les evenements sismiques : 

20 localisation, magnitude, heure, jour de la semaine) et 
des parametres de plus bas niveau extraits des signaux 
par des procedures de type traitement du signal ; 

- un systeme expert flou, dit SEF, prenant 
une decision de fagon independante a partir des memes 

25 informations, et capable d f expliquer sa decision a 
1' utilisateur par 1' intermediaire de regies triees par 
ordre d' applicability a 1' evenement traite ; 

- un reseau neuronal a connexions locales et 
poids partages, dit TDNN, constitue de bancs de filtres 

30 adaptatifs non lineaires, extrayant lui-meme des 
informations pertinentes de representations temps- 
frequence des signaux correspondant a 1' evenement. 

Ces trois branches se configurent 
automat iquement par apprentissage statistique sur des 

35 bases de donnees d' evenements sismiques. 
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• Apprentlssage a partir dL' exewples 

1/ apprentissage a partir d' exemples consiste 
a construire un modele du systeme de prise de decision 
5 par ajustement progressif de parametres a partir de 
donnees . Ce modele doit etre capable d'associer la 
bonne decision (sortie) a un ensemble de donnees 
decrivant le cas trait6 (entrees) . Cela se fait 
progressivement, par presentation iterative des cas 
10 disponibles dans la base d' exemples en entree du 
systeme. Une telle procedure est illustre sur 
1' organigramme de la figure 5. 

Dans 1' invention, le modele qui apprend peut 
etre soit un reseau de neurones artificiels, soit un 
15 systeme expert flou. 

Une fois que le systeme a termine sa phase 
d' apprentissage, ses parametres internes sont figes et 
le systeme est pret a etre utilise. 

20 • Reseau de neurones de type « Perception maltlcouche » 

Un tel reseau de neurones artificiels de type 
« Perceptron multicouche » est un modele particulier de 
reseau de neurones pouvant etre utilise comme systeme 
de prise de decision. II est constitue d' un reseau 
25 d' automates de calcul simples, les " neurones 
artificiels 

Un neurone coirane illustre sur la figure 

6, est une entite constitute d'un vecteur-poids 
Wy={w & -} et d'une fonction de transfert non lineaire <|>. 
30 II admet en entree un vecteur X = {x / } et effectue une 



transformation de ces entrees de type JVy=^X w y x ' 

Par similitude avec le vocabulaire employe en 
neurophysiologie, on dit que chaque entree x± est 
reliee au neurone N J par une connexion synaptique. Un 
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poids synaptique w±j module l'efficacite de cette 
connexion, 

Dans un reseau de neurones artificiels, comme 
illustre sur la figure 7, les neurones sont assembles 
5 en couches successives. One couche est definie comme un 
ensemble de neurones n'ayant pas de connexions entre 
eux, mais pouvant avoir des connexions avec des 
neurones des couches precedentes (entrees) ou suivantes 
(sorties) • En general, on ne connecte que des neurones 
10 de couches successives. 

1/ apprentissage consiste a modifier 

progressivement les valeurs des poids wij jusqu'a ce que 
les sorties du reseau, qui est constitue d'un certain 
nombre de couches de neurones, correspondent aux 
15 sorties desirees. 

Pour cela, on definit une erreur de 
classification que l'on cherche a minimiser. L' erreur 
la plus couramment utilisee est 1' erreur quadratique, 

^SORTIES f . , V2 

definie par E= £ \ z k~ z k aesiree ) . La methode consiste a 

A:=l 

20 faire une descente de gradient sur les poids par 

dE 

1' equation Atv i7 =-a avec a > 0 . Cette equation, 

J 9wij 

lorsqu'on la developpe, donne la formule de correction 
pour chaque poids du reseau. 

25 • Systems expert flou 

Un systeme expert flou est un autre modele de 
systeme de prise de decision. II presente l'avantage 
sur le reseau de neurones de donner une forme 
d' explication de ses decisions. II est constitue d'un 
30 ensemble d' unites de calcul, les M regies d' inference 
floue ". 

Une regie d' inference floue est une entite de 
la forme >x si < premisse > alors < conclusion > " . La 
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premisse est la partie que 1'on cherche a mettre en 
correspondence avec les donnees d' entrees. 

En effet, un ensemble flou est un ensemble 
dont les fronti^res sont progressives, au contraire 
5 d'un ensemble classique, dont les frontieres sont 
nettes. Un element fait done plus ou moins partie de 
chaque ensemble flou. Lorsque les donnees sont de 
dimension 1, un ensemble classique peut etre represents 
par un rectangle (appartenance = 1 a l'interieur, 0 au- 

10 dehors), alors qu'un ensemble flou peut etre un 
triangle, un trapeze, une gaussienne... 

De meme que precedemment , 1' apprentissage 
consiste a modifier progressivement les valeurs des 
parametres jusqu'a ce que les sorties du systeme expert 

15 flou correspondent aux sorties desirees. 

Quatre types de parametres sont calcules par 
apprentissage : la position et la largeur des ensembles 
flous des premisses, les poids des regies et le degre 
deactivation de chaque classe pour chaque regie. 

20 En phase d' utilisation operationnelle, le 

systeme expert flou fournit, en plus de la classe 
attribute a l'evenement sismique, la liste des regies 
applicables par ordre decroissant de pertinence. 
Certaines de ces regies peuvent etre en contradiction 

25 avec les autres, ce qui permet d' examiner les solutions 
alternatives, mais e'est l'agregation du resultat de 
toutes les regies qui fournit le resultat global. 

L' utilisateur a done a sa disposition : 

- la decision brute (seisme, explosion, coup 
30 de terrain) ; 

- la liste des regies applicables ; 

- la liste des regies apportant une 
contradiction a cette decision ; 



B 13340.3 DB 



18 



- la raison de la decision de chaque regie 
(par examen de la coherence entre les donnees et les 
ensembles flous correspondants) . 

Un exemple de regie de decision trouvee par 
5 le systeme est le suivant : 

si (Heure est milieu jie apres-midi) 

et (Latitude est tres j>rochejde 43.5° N) 

et (Longitude est tris j>rochejde 5.5° E) 

et (Magnitude est environ 2.7) 

10 et (Date est de preference samedi) 

alors (avec niveau de confiance = 0,8) 

(s6isme est improbable) 
(explosion est probable) 
(coup de terrain est improbable) 

15 

Dans 1' invention, pour des raisons de 
convergence difficile, on parametre cette descente de 
gradient en introduisant des variables intermediaires . 
Si 1'on veut effectuer une descente de gradient sur un 

20 parametre p avec p=<|>(s), <|> etant une fonction 
derivable, strictement monotone, independante de p et 
des valeurs des exemples servant a 1' apprentissage , on 
a les memes solutions finales en effectuant une 
descente de gradient sur s. L' avantage d' un tel 

25 changement de variable est qu' il devient possible de 
changer la fa<?on d'atteindre la solution, et en 
particulier de faciliter la convergence dans les cas 
difficiles. 

Dans l f invention, on optimise les parametres 

30 suivants : 

(1) La position y des ensembles flous des 
premisses : lorsque la descente de gradient est 
appliquee directement sur ce parametre, on obtient 
generalement une convergence difficile. Ceci s'explique 
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par le fait que la variation de la position y des 
ensembles flous des premisses n'est pas une fonction 
croissante de la distance a l'exemple. On corrige ce 
phenomene en posant x=y/a. 
5 (2) La largeur a des ensembles flous des 

premisses : lorsque les donnees sont structurees en 
groupes de tailles tres differentes, l'algorithme ne 
peut pas converger. En etudiant la variation relative 
Aa/cr, on decouvre qu'elle n'est pas bornee (c'est-a- 

10 dire que rien ne l'empeche de tendre vers des valeurs 
infinies) . Lorsque les donnees sont tres groupees, 
cette variation prend ef f ectivement des valeurs tres 
elevees. Pour avoir une modification relative plus 
faible lorsque les donnees sont rapprochees, on pose 

15 s=ln(2a 2 ). 

(3) Les poids p des regies : c'est le 
parametre le plus difficile a regler. Avec une descente 
de gradient directe, les poids les plus faibles 
diminuent et deviennent negatifs, ce qui leur fait 

20 perdre toute signification et fait diverger 
l'algorithme. On choisit done une fonction d' activation 
positive en imposant une contrainte supplementaire : 
pour differents exemples avec le meme niveau activation 
de la regie, la variation de ce niveau doit etre la 

25 meme si les conclusions sont egales . La consequence est 
que la variation relative des poids des regies doit 
etre constante lorsque les exemples ont le meme degre 
d' appartenance aux ensembles flous. Ceci est realise en 
posant r=ln (p) . 

30 (4) Le degre d' activation d de chaque classe 

pour chaque regie. 

La descente de gradient est done effectuee 
non pas sur y, a et p, mais sur : 
■ x=y/a 
35 ■ s=ln(2cr 2 ) 
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■ r=ln(p) 

Pour d, on n'effectue pas de changement de variable. 
Ces changements de variables assurent une tres bonne 
qualite de convergence et permettent 1'obtention de 
5 systemes-experts flous tres efficaces. 

• Base d r examples efc validation 

La base d' exemples iitilisee doit verifier 
deux principes fondamentaux : 

10 - etre qualitativement representative du 

probleme reel (distribution des exemples conformes a la 
distribution reelle) ; 

- etre cjuantitativement representative du 
probleme (nombre d' exemples suffisant pour constituer 

15 un echantillonnage satisf aisant ) . 

II existe plusieurs methodologies 

d' apprentissage et de validation. Dans la procedure la 
plus simple, on divise la base d' exemples en deux bases 
disjointes : la base d' apprentissage et la base de 

20 test. On entraine le systeme par apprentissage sur la 
premiere et on verifie son bon f onctionnement sur la 
seconde. Une base d' exemples ne verifiant pas 
conjointement les deux proprietes enoncees ci-dessus 
fait courir le risque de conduire a un systeme 

25 incapable de generaliser correctement , c'est-a-dire de 
fonctionner sur des exemples nouveaux, non presentes au 
cours de la phase d' apprentissage . 

Que ce soit pour les reseaux de neurones 
artificiels, pour les systemes-experts flous ou plus 

30 generalement pour tout systeme congu par apprentissage 
statistique sur des donnees experimentales, il est 
primordial d'utiliser une base d' exemples suffisante en 
qualite et en quantite, et de valider le systeme mis au 
point par des procedures rigoureuses . 

35 
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Description detaillee 

• Description d'un evenement sismigue 

Un evenement sismique a identifier peut etre 
decrit par : 

5 - la donnee des signaux en provenance du 

reseau de stations sismiques, ou 

- des proprietes de plus haut niveau, 
directement mesurables ou calculees par des modeles 
geophysiques. On peut utiliser par exemple la 

10 localisation de 1' evenement (latitude et longitude) , sa 
magnitude et le moment ou il est survenu (heure et jour 
de la semaine) . Par exemple les informations de haut 
niveau sont les suivantes : le jeudi 7 avril 1966 a 12 
heures, un seisme de magnitude 1,4 s'est produit par 

15 02°35'06" de longitude est et 49°12'25" de latitude 
nord. 

• Le classifies: neuro-flou 

Le classif ieur neuro-flou (expert CNF) , comme 

20 illustre sur la figure 8, est constitue d'un codage 
neuro-flou des donnees suivi d'un perceptron 
multicouche. II s' applique aux donnees de haut niveau. 

Le codage neuro-flou consiste a associer 
plusieurs cellules de codage a chaque variable (ou 

25 groupe de variables) d' entree, chaque cellule ayant une 
region d' influence modelisee par une f onction 
definissant son mecanisme d' activation . La presentation 
d'un vecteur de valeurs genere alors un schema 
d' activation des cellules de codage associees. 

30 La figure 9 illustre ce mecanisme 

d' activation des cellules de codage. La presentation 
d'une valeur genere un schema d' activation 
correspondant a la reponse impulsionnelle de chaque 
f onction d' activation a la valeur presentee. Les 

35 niveaux de gris affectes aux centres des cellules 
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indiquent leur niveau d' activation compris entre 0 et 1 
(noir : 1, blanc : 0) . 

La figure 10 montre des exemples de schemas 
deactivation generes par la presentation de valeurs 
5 typiques. II s'agit d' un codage de type curseur. La 
valeurs faibles ( respect ivement elevees) activent 
pref erentiellement les cellules de gauche 

(respect ivement droite) . 

L'interet de ce codage est multiple : 
10 - par sa nature meme, il permet de 

representer des donnees incompletes , imprecises ou 
incertaines et de les utiliser ef f icacement pour la 
prise de decision ; 

- par ses proprietes de conditionnement non 
15 lineaire des donnees, il facilite les traitements 
ulterieurs (ici, la classification) . 

Ce codage neuro-flou est realise en plusieurs 
etapes successives : definition des sous-groupes de 
caracteristiques, choix et placement des cellules de 
20 codage affectees a chaque groupe, determination des 
parametres de la region d' influence de chaque cellule . 
Les details de cette procedure sont exposes dans le 
document reference [8] . 

Une fois les donnees codees, elles sont 
25 analysees par le perceptron multicouche qui en calcule 
la classe. 

• Le systems expert flou 

Dans un mode de realisation le systeme de 
30 1' invention comporte une seule branche de traitement 
basee sur un tel systeme expert flou 

Le systeme expert flou (expert SEF) 
s' applique egalement aux donnees de haut niveau. 

Dans la figure 11, on represente un systeme 
35 expert flou a cinq regies (une regie par ligne) . 
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Pour chaque ligne, les cinq colonnes de 
gauche representent les premisses et correspondent a 
cinq entrees : l'heure, la latitude, la longitude, la 
magnitude et la date. Les premisses sont composees 
5 d' ensembles flous en forme de gaussiennes qui 
recouvrent les domaines de la variable d' entree 
conduisant a un renf orcement de l'activite de la regie. 

Quatre types de parametres sont calcules par 
apprentissage : la position et la largeur des ensembles 
10 flous des premisses (colonnes 1 a 5) , les poids des 
regies (colonne 6) , qui permet de specifier le degre 
d' importance de chaque regie dans le processus de 
decision et le degre d' activation de chaque classe 
(seisme naturel, explosion ou coup de terrain) pour 
15 chaque regie (colonnes 7 a 9) . 

A chaque presentation d'un exemple a 
classifier, on calcule le contenu de la colonne 10 et 
de la colonne 6 : 

■ La colonne 10 donne 1' activation de chaque 
20 regie (et permet done d'estimer son adequation avec le 

cas traite) . 

■ La ligne 6 est la synthese des decisions 
des cinq regies et donne la reponse globale du systeme 
expert flou (ici, la decision est « explosion ») . cette 

25 synthese est faite en calculant le barycentre des 
decisions de toutes les regies (colonnes 7 a 9) 
ponderees par le niveau deactivation correspondant 
(colonne 10) . Sur la figure 11, la position du 
barycentre pour chaque classe est symbolisee par un 

30 trait vertical ligne 6, colonnes 7 a 9. 

L' apprentissage est realise en deux temps : 
- une premiere phase consiste a positionner 
les ensembles flous (centres et largeurs), par exemple 
au moyen d'un algorithme dit des C-moyennes floues, 

35 comme decrit dans le document reference [9] ; 
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- une deuxieme phase consiste a realiser une 
descente de gradient sur les quatre types de 
parametres. 

5 • Le reseau neuronal a connexions locales et poids 
paj?ta.ges 

Au contraire des deux branches precedentes, 
le reseau de neurones a connexion locale et poids 
partages (expert TDNN) admet en entree les signaux 

10 sismiques eux-memes et apprend a extraire par 
apprentissage non seulement la procedure de decision, 
mais aussi les parametres discriminants qui vont servir 
de base a cette decision. Ce reseau neuronal est de 
type perceptron multicouche a connexions locales 1 et 

15 poids partages prenant en entr6e les spectrogr amines 
pretraites des signaux sismiques, comme decrit dans le 
document reference [10] . Ces spectrogrammes sont 
obtenus par 1' application d'une transformee de Fourier 
a fenetre glissante sur le signal. 

20 Les figures 12A a 12D illustrent les 

pretraitements successifs appliques a chaque signal 
sismique aboutissant a un spectrogramme final a 15 
bandes f requentielles : la figure 12A illustre le 
signal initial ; la figure 12B illustre le 

25 spectrogramme deduit du signal avec 50 bandes 
f requentielles ; la figure 12C illustre le 
spectrogramme « debruite » ; la figure 12D illustre le 
spectrogramme apres reduction de 50 a 15 bandes 
f requentielles . 

30 Le spectrogramme obtenu est ensuite pretraite 

puis presente en entree d'un reseau de neurones de type 
TDNN. Chaque reseau est specialise dans le traitement 
des signaux enregistres par une station donnee . 

La figure 13 illustre 1' architecture d'un 

35 reseau TDNN specialise dans la classification des 
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spectrogrammes deduits des signaux enregistres par un 
sismometre donne, ce reseau comportant quatre couches 
de neurones. La couche d' entree a des connexions 
locales et des poids partages (4 trames avec un 
5 decalage de 2 trames) avec la premiere couche cachee. 
Celle-ci a de meme des connexions locales a poids 
partages (9 trames avec un decalage de 5 trames) avec 
la seconde couche cachee, totalement connectee a la 
derniere couche. 

10 Les poids partages rendent 1' architecture 

plus robuste aux petits decalages des pointages des 
phases ou de trames manquantes ou erronees. Cependant, 
a cause de la vitesse de propagation des ondes P 
(compression) et S (cisaillement ) , la duree entre 

15 l'arrivee de la phase P et de la phase S varie en 
fonction de la distance entre la station 
d' enregistrement et 1' epicentre de l'evenement, ce qui 
complique 1' apprentissage . La solution adoptee consiste 
a aligner le pointage de la phase P sur la 10 6me trame 

20 et celui de la phase S sur la 60* me trame. 

• La. prise de decision finale 

Pour la prise de decision finale, on suppose 
que toutes les sorties sont comprises dans 1'intervalle 

25 reel [-1,1]. Cette prise de decision consiste en une 
association des reponses fournies par les trois 
branches afin d'augmenter la fiabilite. Elle peut etre 
realisee par un calcul de moyennes arithmetiques sur 
les sorties homologues de chacune des trois branches. 

30 Pour chacune des trois sorties S± du systeme global, on 
a alors : 

..3 i j=\ 

La certitude sur la reponse est evaluee par 
un coefficient calculable seulement si le systeme est 
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en situation de prise de decision (c' est-a-dire s'il y 
a une et une seule sortie strictement positive) . Ce 
coefficient est alors egal a la moyenne des valeurs 
absolues des sorties : 

J 1=1 

K<0.2 : doute 

tfe]o.2,0.4] : prudence 
Ke]0A,0.6] : certitude raisonnable 
ATe]0.6,0.8] : grande certitude 
10 • tf>0.8 : certitude quasi absolue 



Ainsi, on peut avoir par exemple : 



Systfeme 


Classe choisie 


Detail des r£ponses par classe 


Degr6 de certitude , 


1 


Classe 3 


(-0,9 -0,4 +0,8) 


0,7 : grande certitude 


2 


Indetermination 1 ou 3 


(+0,1 -0,5 +0,3) 


Incertitude complete 


3 


Classe 1 


(+0,2 -0,6 -0,2) 


0,3 : prudence 


Fusion 


Classe 3 


(-0,2 -0,5 +0,3) 


0,3 : prudence 



15 

Exemple de mise en oeuvre de 1' invention : 
discrimination d' evenements sismiques reqionaux 
• Localisation de 1' evenement 

La discrimination « evenement naturel / 

20 evenement artificiel » est une etape majeure de la 
surveillance sismique, r6alis6e rapidement des la 
lecture des signaux lors de la phase de depouillement , 
puis affinee au cours de chaque nouveau traitement. On 
evalue a environ sept ans le temps necessaire aux 

25 analystes pour etre veritablement operationnels . Ces 
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analystes etant de veritables experts, il est difficile 
d'exprimer clairement leur methode de raisonnement 
basee a la fois sur des connaissances expertes et sur 
un raisonnement par cas. 
5 La localisation d'un evenement sismique est 

obtenue apres une succession d'aller-retours entre deux 
phases principales : 

- le pointage des differentes phases 
sismiques effectue sur les signaux enregistres par les 

10 stations ayant detecte 1' evenement et le calcul de la 
magnitude ; 

- la localisation elle-meme, realisee a 
l'aide d'un modele mathematique con?u par les 
sismologues . 

15 A 1' issue de la phase de pointage, la 

localisation de 1' evenement peut-etre effectuee a 
l'aide d'un logiciel de simulation connu de l'homme de 
metier. II fait appel a des modeles sismologiques 
contenant des informations sur la vitesse des ondes, 

20 les differents types d' ondes et leur mode de 
propagation, des corrections d' altitude selon les 
stations... Plusieurs hypotheses de localisation sont 
proposees, associees a un degre de coherence des 
donnees. Si celles-ci ne conviennent pas a 1' expert, il 

25 modifie ses pointes de phases puis relance une 
localisation. Ce cycle est reitere jusqu'a l'obtention 
d'un resultat juge satisf aisant . 

La qualite de la localisation depend du 
nombre et de la qualite des stations ayant servi a la 

30 localisation de 1' evenement, ainsi que de leur 
repartition azimutale. Les evenements situes en France 
sont done, en regie generale, mieux localises que les 
evenements etrangers. Pour af finer la localisation dans 
ce dernier cas, les experts accedent regulierement a 

35 des donnees etrangeres . Avec les seules donnees 
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frangaises, la precision de localisation des evenements 
survenant en France est en moyenne de cinq kilometres. 
Elle est estimee a environ un kilometre dans le 
meilleur des cas. 
5 Le tableau II en fin de description donne un 

exemple de resultats fournis par la procedure de 
localisation. La partie superieure resume les resultats 
obtenus : a chaque valeur d'heure origine, de magnitude 
et de localisation (latitude-longitude) est associee le 

10 niveau d' imprecision estime. La partie inferieure 
indique les stations ayant servi £ la localisation et 
les erreurs moyennes quadrat iques (RMS) obtenues en 
fonction de 1'hypothese de profondeur de 1' epicentre. 
Dans le cas des coups de terrain, la profondeur est 

15 arbitrairement fixee a un kilometre. 

• Caracterlsatlon de 1 'evenement 

On s'interesse ici aux evenements sismiques 
locaux et regionaux, c'est-a-dire survenus en France 
20 metropolitaine ou dans les regions limitrophes. Ces 
evenements sont souvent qualifies de proches par 
opposition aux teleseismes lointains (epicentre situe a 
plusieurs milliers de kilometres du capteur) . 

Trois types d' evenements sismiques doivent 
25 etre discrimines : 

- les seismes, evenements sismiques d' origine 
naturelle ; 

- les explosions terrestres (tirs de mines, 
de carrieres, de chantiers...) ou en mer (deminage, 

30 essais d'armes...) ; 

- les coups de terrains correspondant a 
1' ef fondrement d'une couche miniere et lies a 
l r exploitation de la mine. 

L' analyse de l'etat de l'art a montre 1'echec 
35 des approches basees sur une discrimination a partir de 
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signaux . sismiques uniquement. On exploite done 
1' ensemble des donnees disponibles en adoptant une 
approche basee sur la fusion multi-experts et multi- 
sources. La conception du systeme de discrimination 
5 automatique d' evenements sismiques est basee sur trois 
modules : 

- les deux premiers (experts CNF et SEF) 
sont des modules effectuant une discrimination a partir 
des seules donnees de haut niveau, deduites par le 

10 module inverse du laboratoire LDG. Aucun signal 
sismique n' est done directement pris en compte a ce 
niveau ; 

- le troisieme (expert TDNN) est base sur 
1' analyse des signaux sismiques. 

15 

• Donnees utill&ees 

Repartition geographique des evenements 

Les evenements sismiques a analyser se 

repartissent sur 1' ensemble du territoire frangais et 
20 dans un perimetre proche . Les epicentres des evenements 

enregistres par le laboratoire LDG de 1962 a 1996 sont 

illustres sur la figure 14. 




Donnees de haut niveau 

25 Chaque evenement sismique est caracterise 

par les informations suivantes : la date et l'heure 
origine de 1' evenement , la latitude et la longitude de 
1' epicentre, et sa magnitude. 

L'heure et la date sont retenues en raison 

30 des regies regissant les tirs de carrieres ou de mines 
en France, interdisant les tirs nocturnes ou durant les 
wee k- enc ls et jours feries. Neanmoins, des derogations 
sont accordees a certains chantiers, par exemple pour 
eviter les nuisances a la circulation. 
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La magnitude est retenue car les coups cle 
terrain ont une magnitude typique d'apres les 
sismologues (environ 3) . De plus, seuls les seismes 
peuvent fournir les plus fortes magnitudes. Plusieurs 
5 valeurs de magnitude sont prises en compte lorsqu'elles 
sont disponibles. 

La localisation de 1' epicentre, caracterisee 
par sa latitude et sa longitude, est egalement une 
caracteristique majeure. II existe neanmoins certaines 
10 mines situees dans des regions a forte sismicite et 
susceptibles d' occasionner des coups de terrain. 

Donnees de bas niveau 
Les donnees de bas niveau sont les signaux en 

15 provenance des 42 stations sismiques du laboratoire LDG 
(voir la figure 1). Les pretraitements portent 
essentiellement sur la creation des spectrogrammes des 
signaux sismiques, qui sont non stationnaires . Ces 
spectrogrammes sont obtenus par 1' application d'une 

20 transformee de Fourier a fenetre glissante sur le 
signal. Dans un premier temps, le signal, echantillonne 
a 50 Hz, est segmente en trames de deux secondes 
decalees d'une seconde par une fenetre de Hamming. On 
calcule ensuite la densite spectrale d' energie sur 50 

25 bandes f requentielles avec elimination de la composante 
continue. Puis, on applique une transformation 
logarithmique avec soustraction de bruit, dont le 
modele est suppose logarithmique, dans chaque 
bande selon 1' equation mex(ln(l + x)- Meruit)- cr(bruit\o) , ou 

30 ij (bruit) et a(bruit) correspondent a la moyenne et a 
l'ecart-type du bruit estimes sur une periode 
anterieure au pointage de l'onde P. Enfin, on reduit le 
nombre de bandes f requentielles de 50 a 15 par une 
compression pseudo-logarithmique des hautes frequences. 

35 
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• Resultats obtenus 

Le systeme decrit ci-dessus classifie les 
5 6venements sismiques regionaux frangais avec les 
performances suivantes : 

- 86 % pour les seismes ; 

- 91 % pour les explosions ; 

- plus de 99% pour le coups de terrain. 
10 Les performances globales sont d' environ 90%. 
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TABLEAU II 



nombre de stations utilisees 11 - Res dans bulletin: oui 
preloc. lat: 43.2 Ion: 5.4 ho: 19 19 53,1 



region: 11 to S Aii-en-Provence (13) 

• c de toit 1 

heme origine: 15 19 53.8 +/- .3 
latitude: 43.43 ♦/- 1.7 

longitude: 5.43 ♦/- 2.5 
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RE VEN D I CAT IONS 



1. Systeme d' intelligence artif icielle de 
classification d'evenements, d'objets ou de situations 

5 a partir de signaux et de parametres discriminants 
issus de modeles, caracterise en ce qu' il comprend au 
moins une branche de traitement comprenant un systeme 
expert flou prenant une decision a partir de proprietes 
de haut niveau et de parametres discriminants de plus 
10 bas niveau extraits des signaux par des procedures de 
type traitement de signal, et capable d' expliquer sa 
decision a 1' utilisateur par 1' intermediaire de regies 
triees par ordre d' applicability. 

2. Systeme selon la revendication 1, dans 
15 lequel, dans le systeme expert flou on effectue une 

descente de gradient sur les parametres : 
x = y/o 
. s = In (2a 2 ) 
- r = ln(p) 
20 - d 

avec : 

y : position des ensembles flous des 
premisses 

a : largeur des ensembles flous des 
25 premisses 

p : poids des regies 
. d : degre deactivation de chaque classe 
pour chaque regie 

3. Systeme selon la revendication 1, qui est 
30 un systeme multi-experts constitue d'au moins deux 

branches de traitement independantes, se configurant 
automatiquement par apprentissage statistique sur des 
bases de donnees, ayant des proprietes particulieres et 
fusionnees par un systeme decisionnel de haut niveau. 
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4. Systeme selon la revendication 3, dans 
lequel une branche comprend un classifieur neuro-flou 
prenant ses decisions a partir de propriety de haut 
niveau et de parametres discriminants de plus bas 

5 niveau extraits des signaux par des procedures de type 
traitement de signal. 

5. Systeme selon l'une quelconque des 
revendications 3 ou 4, dans lequel une branche comprend 
un reseau neuronal a connexions locales et poids 

10 partages constitues de bancs de filtres adaptatifs non 
lineaires, extrayant lui-meme des informations 
discriminantes de representation temps-f requence des 
signaux correspondants . 

6. Systeme selon l'une quelconque des 
15 revendications precedentes, qui est un systeme de 

classification d' evenements geophysiques . 

7. Systeme selon la revendication 6, dans 
lequel les proprietes de haut niveau sont la 
localisation, la magnitude, l'heure et la date. 
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